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Abstract 

The flat-space limit of the one-loop effective potential for SO(IO) GUT the- 
ories in spatially flat Friedmann-Robertson- Walker cosmologies is applied to 
study the dynamics of the early universe. The numerical integration of the 
corresponding field equations shows that, for such grand unified theories, a 
sufficiently long infiationary stage is achieved for suitable choices of the initial 
conditions. However, a severe fine tuning of these initial conditions is neces- 
sary to obtain a large e-fold number. In the direction with residual symmetry 
SU{4:)ps SU{2)l SU{2)ji, one eventually finds parametric resonance for 
suitable choices of the free parameters of the classical potential. This phe- 
nomenon leads in turn to the end of infiation. 
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I. INTRODUCTION 



Semiclassical effects in quantum field tfieory are at tfie very lieart of many exciting de- 
velopments in modern theoretical physics, e.g., radiative corrections to the Casimir force 
[1], trace anomalies [2], one-loop effective action [3,4], symmetry breaking in the early uni- 
verse [5-7]. In the latter class of phenomena, the one- loop effective potential for grand 
unified theories (GUT) in curved backgrounds is used in the field equations to determine 
the symmetry-breaking pattern which is relevant for the familiar low-energy phenomenology 
of particle physics [6,7]. In particular, in a previous paper by the first two authors, the 
one-loop effective potential for SO (10) GUT theories in de Sitter space was derived for the 
first time [8]. The main motivation of our work was the analysis of spontaneous symmetry 
breaking in the early universe when GUT theories in good agreement with low-energy phe- 
nomenology are studied. Unlike SU{5) theories, SO{10) models make it possible to allocate 
a whole fermionic family in the single irreducible spinorial representation, and predict longer 
nucleon lifetimes in agreement with experimental data. This happens because lepto-quarks 
responsible for nucleon decays have higher masses (with respect to SU{5)), by virtue of an 
intermediate symmetry group. The latter lies in between the highest scale Mx = 10^^ GeV, 
and the electroweak scale. Further details can be found in Refs. [8,9]. 

On considering the particular mass matrix relevant for the SU{3)c ® SU{2)i ® U{1)y 
symmetry-breaking direction, our main result was the numerical proof that, even when 
curvature effects are no longer negligible, the early universe can only reach the SU (4) ps ® 
SU{2)l® SU{2)r absolute minimum [8]. 

A naturally occurring question is whether the analytic formulas obtained in Ref. [8] 
can be used to gain a better understanding of other phenomena occurring in the early 
universe. Indeed, it is well-known that, according to the cosmological standard model based 
on Einstein's general relativity, the early universe is spatially homogeneous and isotropic 
[10,11], and gravity couples to the energy-momentum tensor of matter. Under the previous 
assumptions of symmetry and homogeneity (cf. Ref. [12]), the space-time metric can be 
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locally cast in the form [10] 



g = -dt 0dt + a?{t) dx0dx + f{x) {d6i d0i + sin^ Oi de2 dO-^ 



(1.1) 



where a{t) is the cosmic scale factor, 6i e]0,7r[, 62 e]0,27r[ and, denoting by k the constant 
curvature of the three-dimensional spatial sections, /(x) = sinx, X or sinhx if /c = 1, 0, — 1 
respectively. Here we focus on spatially flat Friedmann-Robertson- Walker (FRW) cosmolo- 
gies, for which k — Q. 

The aim of the present work has been therefore to study the relevance of SO(IO) GUT 
models for gravitational physics within the more general framework of FRW cosmologies. 
In this case, the analysis performed in Ref. [13] enables one to use, at least for typical 
GUT-inflationary models characterized by a phase-transition temperature < 10^^ GeV , the 
flat-space limit of the Coleman- Weinberg one-loop effective potential [5]. What happens 
is that the one-loop field equations in FRW cosmologies are an involved set of integro- 
differential equations [13]. In scalar electrodynamics, the nonlocal correction to the Coleman- 
Weinberg effective potential results from the coupling of a scalar field 0c to the gravitational 
background and may lead to dissipative effects in the early universe. More precisely, nonlocal 
terms in the field equations may be approximated very well by a function of the form [13] 



where ^4 is a positive-definite function of the conformal time f. The first term in round 
brackets leads to a further quantum correction to the energy density of the scalar field (for 
slowly varying A), while the second term is purely dissipative [13]. However, nonlocal terms 
do not modify the pattern of minima found in static de Sitter space. Thus, for numerical 
purposes, even in the nonAbelian case, curvature effects on the one-loop potential are very 
small, and can be neglected at the GUT energy scale. This is the basic approximation on 
which the present paper relies. 

In Sec. II the scalar dynamics driving the inflationary phase is described. Section III, 
relying on Ref. [8], studies SO(IO) GUT models with Higgs field in the 210-dimensional 
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irreducible representation. In the SO{10) internal space, spherical coordinates for the Higgs 
field are used to obtain a convenient parametrization of the classical part of the potential. 
The flat-space limit of the one-loop effective potential is then studied in Sec. IV. Semiclas- 
sical field equations are studied in Sec. V, and their numerical solution is found in Sec. VI. 
The reheating process is then briefiy described in Sec. VII, while the concluding remarks 
are presented in Sec. VIII. 



II. INFLATION FROM SCALAR-FIELD DYNAMICS 

In the infiationary models driven by scalar fields, the universe contains different species 
of matter, i.e., radiation and scalar particles (inflatons). The latter, characterized by a 
higher energy level with respect to the others, completely determine the dynamics of the 
universe, making all remaining contributions negligible. A theoretical scheme in which 
fundamental scalar particles naturally occur (Higgs fleld) and lead to symmetry breaking are 
the GUT theories. In this case the energies involved are compatible with the levels required 
for an efficient inflation. Thus, we consider an inflationary model in which the inflaton 
coincides with the Higgs fleld of a GUT theory, corresponding to the highest energy scale 
of spontaneous symmetry breaking. In particular, we choose the most relevant ones which 
are the SO(IO) GUT theories with Higgs fleld belonging to the 210-dimensional irreducible 
representation ( 210 ) [8]. 

Within this framework (we denote hereafter by the flelds corresponding to scalar parti- 
cles), and assuming homogeneous fleld conflgurations = (f>{t), the semic/assica/ Lagrangian 
reads 



C = 6d'aF{(j)) + 6da'F'{(j))(j) + a 



(2.1) 



where the dot stands for d/dt, the prime for 5/d4> and F{4>) denotes the arbitrary couphng 

r ■■ " 2 "I 

of (f) to the gravitational background with scalar curvature R — 6 f + ^ ■ Eq. (2.1) 
the classical potential has been replaced by the one- loop effective potential V^^\(f)), i.e., the 
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gauge degrees of freedom have been integrated over. Note also that the first two terms in 
(2.1) yield the coupling —a^F{(j))R, after integration by parts in the action. 
Starting from (2.1) one gets the two semiclassical equations of motion 



a ^ 1 /d 
a 2 \a 



+ 



aF' ■ IF'-- IF"-. 



a F 



2 F 



2 F 



1 02 1 ^ 
8^ + 4^ = °' 



(2.2) 



and 
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i^'(0) + 
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= . 



(2.3) 



In general, for a nonminimally coupled real scalar field, F(0) = (l/167rG) + (^/2)02. in 
our case, in the light of the results obtained in Ref. [13], we neglect all contributions of the 
curvature to the scalar potential, and hence we can assume minimally coupled scalar fields 
(for which ^ vanishes). 



III. SO (10) GUT MODELS WITH HIGGS FIELD IN THE 

210 IRREDUCIBLE REPRESENTATION 

In a gauge unifying theory like 5*0(10), the Higgs field (fundamental scalar particles) 
belongs to one or more irreducible representations (hereafter referred to as IRR's) of the 
gauge group, and its dynamics is ruled by a Higgs potential. These particles provide the 
correct residual symmetry for the model in the low-energy limit, through a spontaneous 
symmetry-breaking mechanism. To study the inflation corresponding to the highest energy 
scale (Mx) of spontaneous symmetry breaking (SSB), we can consider the only contribution 
of the IRR responsible for the SSB at that scale. 

In the present case, we consider the most general renormalizable Higgs potential con- 
structed by using a massless and minimally coupled IRR 210 only, which is obtained by the 
completely anti-symmetrized product of four different lO's as 

^abcd = N ii^^^®^^® Pc®ad] , (3.1) 
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where is a normalization constant. The 210 IRR has four independent quartic invariants, 
i.e., 11011^ and three non-trivial invariants, hence the Higgs potential we are going to construct 
is a function of these [8] 

V{<P) = Vo + g, 11(00)4511' + ^2 \m)mf + 9s Il(00)io5oir + A (3.2) 

where gi, g2, g^ and A are arbitrary coefficients, and Vq is an arbitrary constant. We will see 
in due course what is the meaning of this constant and how it can be fixed. Unfortunately, 
in view of the technical difficulties to express (3.2) in terms of the 210 degrees of freedom, we 
restrict our analysis to the only directions invariant under the subgroup SU {3)c®>SU (2)^(8) 
U{1)y- This restriction, however, remains relevant for the aims of this paper, since it 
leads to the correct electroweak phenomenology at low energies for the model. The most 
general singlet 4>o with respect to the group SU (3)c ® SU {2)l®U (l)y contained in the 210 
representation is [8] 

-j= sin 9 sin </? (01278 + 03478 + 05678 + 01290 + 03490 + 0569o) 



+ ^ sin 61 COS LP (01234 + 03456 + 05612) + COS 9 (0789o) , (3.3) 



l|0o|| 76 

1 

73 

where 9 e [0, tt] and ip e [0, 27r[. In particular, by varying 9 and (p in their ranges, one gets 
the following residual-symmetry groups [8]: 

^ = 0,7r and/or ^p ^ Q,t: ^ SU{'i)c ® SU{2)l® SU{2)r®U{1)b-l , (3.4a) 



9 = ^ and ^ = Q,t:^SU{Z)c®SU{2)l®SU{2)r®U{1)b-l^D, (3.4b) 



9^0,Tr^ SU{4)ps ® SU{2)l ® SU{2)r , (3.4c) 
9 — arctan(3) and (p — arctan(-\/2) or (3.4d) 



9^ - arctan(3) + t: and Lp = arctan(\/2) +t: ^ SU{h)®U (1) , (3.4e) 

6 



otherwise ^ SU{3)c ® SU{2)l U{1)t,^ ® U{1)b-l , 



(3.4f) 



where T^ji is the ^-component of the SU{2)ji group. 

Inserting (3.3) into (3.2) one finds the tree-level potential 

S 



T/(0o) = Fo+(|/a + jf, + ^fs + (A - 5)) ll^of , 



where [81 



(3.5) 



a 



945 V 



108^1 + 28^2 + 140^3) , 



(3.6) 



^^35^^' 



(3.7) 



(3.8) 



n 2 3 

fa = sin^ 6 + sin^ 6 sin^ 99 2 sin 6' cos 99 + a/3 cos 6 + - sin^ ^ sin*^ ip , 



(3.9) 



= sin^ 9 



n 2 



COS 9 COS </? H — -j= sin sin^ (p 
v3 



+ sin^ sin^ (/? cos'' (p + fa , 



(3.10) 



2 sin^ 6* cos^ ~ \ ^ ^^"^ - 3 cos^ 9 



+ 30A - 25/a 



(3.11) 



Since in the following analysis 5 is always negative and a may take negative values, the 
tree-level potential (3.5) is unbounded from below, unless we impose the restriction [8] 



max 9 V ' max 



(3.12) 



Note also that contributions proportional to a cubic term in the potential are set to zero, 
since we are assuming — > —0 invariance of our model. 
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IV. ONE-LOOP EFFECTIVE POTENTIAL 



As shown by Coleman and Weinberg [5] , the effects of quantum fluctuations on the scalar 
potential, due to both gauge bosons and self-interactions, may account for the spontaneous 
symmetry breaking occurring in gauge theories without assuming from the beginning the 
presence of negative quadratic terms for the Higgs held. These phenomena are best tackled 
in terms of the one-loop effective potential [5-8,13,14]. This provides, in our case, the 
appropriate tool for studying the phase transition due to scalar dynamics. 

The one-loop effective potential V^^^ corresponding to the classical expression (3.5) in 
the flat-space hmit is obtained in Ref . [8] . On deflning 



hi = cos^ 9 + sin^ 9 - sin^ ^ + - cos^ 9? 




(4.1) 



/i2 = cos^ 9 + sin^ 9 - sin^ 9? + - cos^ (p 




- sin 9 cos 9 sin 



(4.2) 



/i3 = - sin^ 9 sin^ (p , 



(4.3) 



/i4 = - sin^ 9 , 



(4.4) 



(4.5) 



one gets 




3g^ 



\\M'[lhl{3-Hhl)) + lhl(3-Hhl)) 



87r2 




(4.6) 



where Q is the gauge coupling constant and /x is the renormalization mass. 
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V. SEMICLASSICAL FIELD EQUATIONS 



After defining the dimensionless time r = ^t/Q the Higgs field 0o can be seen, in the 
5'O(10)-space, as the position vector in i?^, given in spherical coordinates as = {Q / lj)4'o = 
y iy. Thus, by taking its derivative with respect to r one has 



d ~ dy 



dO 



dip 



(5.1) 



With our notation, it is convenient to introduce a dimensionless expression for (4.6) obtained 
multiplying it by ^^//x^, i.e., 

5 



V^(^) (y, ^, <^) ^ g/, + 2/^ + + ( A - <5)) y' 



sg' 4 



(5.2) 



As far as the constant Vq is concerned, it has to be fixed by requiring that a vanishing 
potential energy should correspond to the absolute minimum for V^^\ By denoting with 
ym, and Lpjn the absolute-minimum coordinates, we have V'^^\ymiOm-i'4^m)—^- Moreover, 
ym makes it possible to determine also bearing in mind that the spontaneous symmetry- 
breaking scale is Mx, which is fixed for the particular model by the low-energy predictions, 
and then — Mx/ym- 

In this scheme, the semiclassical equations for 0o [15] take the form (from now on, the 
dot denotes differentiation with respect to r) 



y = y 



r + sin^^0^ -?,Ky 



dy 



^ = -2^^ + sin ^ COS ^ 0?^ - 3X^./p^ 

y ^ y^ o6 

V y^sm Bo(fi 



y sin 9 
With our notation, K = ^87rG/^V3e?2, and 

1 



(5.3) 
(5.4) 
(5.5) 



P4> 



(5.6) 



Inflation eventually ends when, by virtue of decay mechanisms, the energy stored in the 
scalar configuration is released to hghter degrees of freedom (radiation). This is the so- 
called reheating process, and the knowledge of the value of r, say r/, for which this actually 
happens, enables one to determine the total e-fold number N{Tf) of the inflationary model 
by solving the differential equation 



A brief description of the reheating mechanism for our particular model is presented in Sec. 



A necessary requirement to perform a numerical analysis of the differential equations 
obtained in Sec. V, is a good knowledge of the absolute minimum of the potential (5.2). 
As already stated in Sec. Ill, if the inequality (3.12) is satisfied, the effective potential is 
bounded from below, but unfortunately, this does not ensure that quantum corrections to 
the tree-level potential will make it possible to fix the value of the field at the minimum, 
Um, so that it is sensibly different from zero. Thus, one has to choose among the values of 
a, 7, 5 and A satisfying (3.12), the ones for which a symmetry breaking actually occurs [8]. 

In Tab. 1, we report for two choices of the free parameters of the classical potential (3.5) 
the corresponding values of 9m, fm, Um and hence the values of /i and Vq. The parameters 
a, 7, 5 and A have been chosen in such a way that small couplings are achieved, and the 
minimum of y is of order 1 . 

In Fig. 1, we plot the solution for y (obtained by means of the NAG Library routine 
D02CAF) corresponding to the initial conditions relevant for the SU{4:)ps®SU (2)l®5'C/(2)r 
symmetry- breaking direction (see (3.4c)), i.e. 9{0) = 0,7/(0) = 10~^,^(0) = 0,y{0) = 
4. X 10~^. Note that a vanishing value of ^(0) and ^(0) is the one resulting from our partic- 
ular choice of residual-symmetry direction, along which the inflationary dynamics evolves. 




(5.7) 



VII. 



VI. NUMERICAL ANALYSIS 
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Moreover, a very small value of y(0), which approaches zero, reflects our choice to start 
from a singlet state which has a complete SO{10) symmetry (hence y{0) = 0). Of course, a 
y{0) value which differs from zero is only taken for numerical convenience, but our results 
are essentially independent of y(0), providing this is very small. For numerical purposes, it 
has been thus convenient to re-express the Eqs. (5.3)-(5.5) in terms of the variables defined 
in Eqs. (A2)-(A4) of the Appendix. However, we keep using the y, 9, ip parametrization 
in our paper, since it makes it easier to describe symmetry breaking. What happens is 
that y characterizes the dynamics which is independent of the symmetry-breaking direction, 
whereas and (p define the various symmetry-breaking directions. After a slow-roll phase 
which is not shown in Fig. 1 for the sake of clarity, y starts increasing until it reaches a 
region where it oscillates in the neighboorhood of a relative minimum. The corresponding 
value of 9 remains equal to 0. Interestingly, the time necessary to reach the region of rapid 
oscillation for y is of order 4. x 10~^^ sec, and the corresponding e-fold number (see (5.7)) is 
~ 100. However, such a large e-fold number is only achieved with the help of a severe fine 
tuning of the y(0) value. This is indeed a peculiar property of Coleman- Weinberg potentials, 
which have a relative maximum at y = 0. Further details about the oscillating phase and 
the end of infiation will be given in Sec. VH. 

In Fig. 2, y is plotted against r when the choice of a, 7, 5 and A described on the second 
line of numerical values of Tab. 1 is made. In such a case, the following initial conditions 
are chosen: ^(0) = 0, ?/(0) = 10-^ (^(0) = 7r/4, ^(0) = a/2 x 10-^ ?/(0) = 3. x 10-^ 0(0) = 0. 
In this case, however, y oscillates in a neighboorhood of the absolute minimum (cf. Fig. 
1). The same happens for and whose evolution is not shown for brevity. Moreover, a 
fine tuning of the initial conditions leads again to an e-fold number of order 100, and the 
duration of the inflationary stage is the same as in Fig. 1. 

The cases described in Figs. 1 and 2 are indeed also relevant for the analysis of the 
reheating stage of the early universe, as shown in Sec. VII, where a brief outline of such a 
phenomenon is presented. 
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VII. REHEATING 



When the energy stored in the scalar-field configuration is released to the relativistic 
degrees of freedom (hereafter denotes the corresponding energy density) and begins to 
dominate the total energy of the universe, the evolution undergoes the so-called reheating 
phase. Such a phase, which is the necessary last stage of the infiationary dynamics, is as 
important as the exponential expansion itself, since it ensures the end of the exponential 
growth and the reheating of the universe. In fact, if the radiation energy dominates the total 
energy of the system, through the equation of state, the energy density and the pressure pn 
satisfy the condition pR + 3pij = ApR > 0, which implies a < 0. 

According to the recent models of reheating, this process can be divided into three stages 
[16-18]. In the first stage, the scalar field decays into massive bosons via parametric 
resonance. The second stage consists in the decay of previously produced particles, and 
the last stage leads to thermalization (for our purposes, only the first two stages will be 
considered) . 

As far as the parametric resonance is concerned, it can be studied starting from the 
equation for quantum fluctuations of a scalar fleld uj, quadratically coupled to 4>. As shown 
in Ref. [16], on setting z — m^t, such an equation may be cast in the form 

uj'l-^[A{k)-2qcos{2z)\uJk^Q , (7.1) 

where the prime denotes differentiation with respect to z. Moreover, A{k) — k'^/m^a^ + 
2q {k/a being the physical momentum), and q = (y'^$^/4m^. In the model considered in 
Refs. [16-18], $ is the amplitude of oscillations of the fleld 4>, and gi is a small coupling 
constant. Interestingly, an exponential instability of the solutions exists and it can be 
interpreted as a rich particle production. This phenomenon is best tackled by studying the 
stability/instability chart of the Mathieu equation [16-18]. 

In particular, the fleld cu may be given by the fluctuations of the scalar fleld itself. In 
such a case, one has to start from Eq. (A9) of the appendix, where the effect of 3iJ^ is 
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neglected, following Ref. [16]. Parametric resonance is only achieved if M{zi, 2;°, z^) therein 
is positive for some values of a, 7, 5, A, in the neighboorhood of which the scalar field starts 
oscillating. In Figs. 3 and 4, we plot the function M against 9 and ip, for the values of a, 7, 5 
and A reported on the first and second line of numerical values of Tab. 1, respectively. In 
each figure, the closed curve represents those values of 9 and ip for which M vanishes. Within 
the region bounded by such a curve, M takes negative values, whereas it is positive outside. 

Interestingly, in Fig. 3, 9m and (pm he in the region of parametric resonance, i.e., outside 
the region where M < 0, while the converse holds in Fig. 4. Of course, since Fig. 1 
corresponds to a case when 9 remains equal to zero, parametric resonance is indeed achieved 
because, for this particular choice of parameters, M{y,9 — 0, (p) is positive, as shown in Fig. 
3. 

Thus, our analysis shows that the inflationary dynamics with residual symmetry 
SU{4:)ps (8) SU{2)l (g) SU{2)ji (see (3.4c)) leads to slow roll and also to parametric reso- 
nance. These properties add evidence in favour of SO(IO) GUT models being physically 
relevant, although one should bear in mind that fine tuning problems remain. However, the 
case described by Fig. 4 is not, by itself, ruled out. One has instead to consider a scalar 
field uj which is not given by the fluctuations of cj) itself [16-18]. 

We now briefly consider the second stage of the reheating process. Once that the quantum 
fluctuations of the 210 IRR have been produced via parametric resonance, they have to decay 
into relativistic degrees of freedom, whose mass is negligible with respect to the Mx scale. 
Note that the decay of the 210 representation into a fermion-antifermion pair cannot occur, 
since it makes it necessary to introduce a Yukawa coupling of the 210 to fermions, which 
would lead in turn to an undesirable mass term for fermions of order Mx- On similar 
ground, any gauge boson coupled to 210 would acquire mass at the scale Mx, and hence 
such a couphng cannot lead to decays into relativistic degrees of freedom. One is thus left 
with decays into other Higgs particles of smaller mass. Indeed, the presence of lighter Higgs 
fields in GUT theories is particularly evident in SO (10), where the spontaneous symmetry 
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breaking which leads to the standard electroweak model occurs in two steps at different mass 
scales: i.e., 

SO{10) SU{3)c (8) SU{2)l ® U{1)y , (7.2) 

where Mx is of order 10^^ — 10^^ GeV to be compatible with the lower limit on proton 
decay, G is one of the intermediate symmetry groups appearing in Eqs. (3.4a)-(3.4f), and 
the scale Mr is the one relevant for neutrino physics. Note that every breaking phase 
in (7.2) is mediated by a different scalar field belonging to a IRR of S'O(IO). A reahstic 
model can be constructed for example by considering in addition to the 210 (0), the scalar 
fields (■0, ■0) belonging to the representation 126 126 . and two ten-dimensional irreducible 
representations (x) which mediate the electroweak symmetry breaking [19] 

SO{10) SU{3)c SU{2)l ® U{1)y ^ SU{3)c ® U{1)q . (7.3) 

In this model, to obtain the above symmetry-breaking pattern, one has to replace the simple 
tree-level potential of Eq. (3.2) by a much more complicated expression V{4>,ip,ip,x) (see 
Eqs. (7)-(ll) of Ref. [19]). Still, as far as the inflation at the highest energy scale is 
concerned, only the terms of the potential containing uniquely the 210 representation with 
the associated scalars are relevant. This is why we neglected in (3.2) the contributions 
resulting from the other representations. The other terms, however, are important when the 
second stage of reheating is considered, since they can mediate the decay of the massive Higgs 
in the lighter ijj, ip and %. The interacting terms between different scalar representations are 
the ones actually relevant in that they do not provide mass for light particles, while decay 
processes are allowed. Such interacting terms occur in quartic form, involve the whole set of 
representations of 5*0(10), and hence, being linear in 0, do not provide mass for the lighter 
Higgs particles (Eq. (11) of Ref . [19]). 
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VIII. CONCLUDING REMARKS 



Our paper has studied the relevance for inflationary cosmology of a nonsupersymmet- 
ric GUT theory which is consistent with the available data on nucleon lifetimes. This is 
the SO{10) model in the 210-dimensional irreducible representation (see Refs. [9,20] and 
literature cited therein). 

Starting from a quartic tree-level potential in the case of minimal coupling, we have 
studied the flat-space limit of the semiclassical field equations, where the one-loop effective 
potential contains the logarithmic terms which result from the Coleman- Weinberg method 
for the evaluation of radiative corrections [5]. Following Ref. [8], when the Higgs scalar 
field belongs to the 210-dimensional irreducible representation of SO {10) (this corresponds 
to the highest energy levels), attention has been restricted to the mass matrix relevant for 
the SU{3)c ® SU{2)l ® U{1)y symmetry-breaking direction, to agree with the low-energy 
phenomenology of the standard model of particle physics. 

The main result of our investigation is that the inflationary dynamics with residual sym- 
metry SU{A)ps ® SU{2)l ® SU{2)r leads to a sufficiently long inflationary stage, and then 
a reheating process occurs via parametric resonance and subsequent decay of the particles 
produced previously. Such processes are three-body decays of 210 into the massless compo- 
nents of the 126 . 126 and K) representations, and result from the mutual quartic coupling 
of all irreducible representations. The intermediate symmetry SU (4) ps ® SU (2) l ® SU {2)ji 
has been extensively studied in the literature on particle physics, since it occurs in the most 
promising S'O(IO) models [9,20]. Thus, a natural and deep link between low-energy phe- 
nomenology, grand unification, inflationary cosmology and physical processes in the very 
early universe seems to emerge from our work. 

Of course, we have only studied some particular yalues of the parameters of the tree- level 
potential. Although particle physics and cosmology lead to restrictions on such parameters, 
our choices are by no means exhaustive. It now appears interesting to get a more quan- 
titative understanding of the reheating process outhned in Sec. VII. Moreover, from the 
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point of view of perturbative properties of quantum field theory, it is necessary to study in 
detail the nonlocal contributions to the semiclassical field equations. These arise already in 
scalar electrodynamics [13], by virtue of the coupling of the scalar field to the gravitational 
background, and are receiving careful consideration since they are related to dissipative and 
nondissipative phenomena in the early universe [13,21]. It also appears interesting to study 
the corrections to the effective potential resulting from a suitable generalization to 5*0(10) 
models of the technique described in Ref . [22] . 
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APPENDIX: 



To obtain the equation for the quantum fluctuations of the 210 IRR, we follow the 
notation of Ref. [8], for which the most general singlet 0o with respect to the gauge group 
SU{3)c (8) SU{2)l ® U{1)y is written down in Eq. (3.3). Within this framework, the 
Lagrangian for scalar flelds only, which are assumed to be spatially homogeneuos, reads 



1 / dzi 
2 



1 / dz2 
2 



1 / dz^^ 
2 



+ E 9'^^i,iZi,j - V{zi, Z2, Z3) , (Al) 

1=1 



where there is summation over the repeated indices i and j. We now take the one- loop 
expansion of (Al) by writing Zi = z^ + Szi, and we make the ansatz 5zi = 5z for simplicity. 
The three independent degrees of freedom of the scalar field turn out to be 



Zi = \\4>o\\ sin 6* cos , 
^2 — Il0o|| sin^sin</7 , 
^3 — \\4>o\\ cos 6* . 



(A2) 
(A3) 
(A4) 
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Thus, the tree-level potential (3.5) can be re-written in terms of (A2)-(A4) by pointing out 
that 



m? = {zlf + {zlf + {zl)\ (A5) 

3 



- {z'i4 + + i4z'.r + /alitor , m 

fsUoW = (30/, - 25/„) ll^of + - ^ - ^(4?) ■ (A8) 



One thus gets the following equation for the quantum fluctuation: 



(PSz dSz 
3H- 



dt^ dt 



- + M{zlzlzl) 

(Jj 



6z = , (A9) 



where k = yk'^, M(z^, Z2, z^) is defined as 



M{z^^, zl z^s) = I (ciio + cioi ^?^3° + Con z^^ + C200 {zif + C020 {Z2Y + C002 (2:3)') , 



(AlO) 



and the Cabd coefficients are given by 



{3 + V3) a + 7(^l + ^-f+^{2 + y/3) 5 + 8X, (All) 

^a+(l + ^y+fV3S + 8X, (A12) 

(^ + V3)a+(^3 + -^j7 + f V3 5 + 8A, (A13) 
3 7 40 

C200 = -a + -^ + —5 + lOA , (A14) 



Clio 



ClOl = 



Con 



C020 



+ ^)« + ^(f + ^)7 + f (5 + 2v^)5 + 10A, (A15) 

3 

C002 = -a + 7 + lOA . (A16) 



1 /39 

2 V8 
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Table 1. 



a 


7 


5 


A 






ym 


(10^6 GeV) 




-3. X 10-2 


5. X 10-^ 


-5. X 10-6 


1. X 10-2 


1.31 


0.94 


1.68 


0.34 


2.49 X 10-2 


-1. X 10-^ 


5. X 10-^ 


-5. X 10-^ 


3.39 X 10-5 


1.70 


1.35 


2.09 


0.27 


3.50 X 10-2 



Figure captions: 

FIG. 1. The solution of Eq. (5.3) is plotted against r, for the first set of values of a, 7, 5, A 
shown in Table 1. 

FIG. 2. The solution of Eq. (5.3) is plotted against r, for the second set of values of a, 7, 5, A 
shown in Table 1. 

FIG. 3. The function defined in Eq. (AlO), divided by y^, is plotted against and for 
the first set of values of a, 7, 5, A given in Table 1. 

FIG. 4. The function defined in Eq. (AlO), divided by y'^, is plotted against 6 and for 
the second set of values of a, 7, 5, A given in Table 1. 
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Fig.l 




Fig. 2 




